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ABSTRACT 

This  report  presents  mathematitvi  techniques  for  calculating  dielectrics 
which  give  rise  to  prescribed  reflection  coefficients  in  certain  problems 
of  one-dimensional  electromagnetic  propagation.  The  techniques  are,  in 
principle  at  least,  exact.  They  are  based  on  the  use  of  an  equation  of 
the  GeTfand-Levitan  type  for  the  one-dimensional  Schrodinger  equation. 
Although  no  practical  applications  are  given,  it  is  hoped  that  this  report 
will  encourage  the  use  of  newer  techniques  in  synthesis  problems. 
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PROPERTIES  OF  DIELECTRICS 
FROM  REFLECTION  COEFFICIENTS  IN  ONE  DIMENSION 


I.  INTRODUCTION 

This  report  is  based  on  a  series  of  lectures  by  the  senior  author  and  has  been  only  slightly 
revised  for  this  written  presentation.  Hence,  the  report  is  more  informally  organized  and  is 
probably  much  longer  than  would  be  the  case  if  a  thorough  revision  had  been  made.  Perhaps, 
however,  the  great  detail  will  enable  a  novice  in  the  inverse  problem  to  see  through  the  math¬ 
ematics  more  clearly. 

The  report  is  concerned  with  the  mathematical  problem  of  determining  the  character  of  a 
dielectric  scatterer  from  the  reflection  coefficient  in  one -dimensional  electromagnetic  scattering. 
Closely  connected  with  this  problem  is  the  "synthesis  problem"  in  which  one  wishes  to  construct 
scatterers  with  prescribed  reflection  or  transmission  characteristics.  It  is  not  generally  known 
that  powerful  techniques  are  available.  The  objective  of  this  report  is  to  present  these  techniques 
in  a  comprehensive  fashion  for  the  benefit  of  those  who  wish  such  information. 

The  technique  which  we  use  is  to  map  tlie  electromagnetic  equations  into  the  one-dimensional 
Schrddinger  equation.  We  then  use  the  Gel 'fand -Levitan  algorithm  to  compute  the  potential  and 
the  wave  function  from  the  scattering  coefficient.  The  results  for  the  electromagnetic  problem 
are  then  obtained  from  the  mapping. 

The  usual  procedure  for  finding  dielectrics  v/ith  prescribed  scattering  properties  is  to  con¬ 
sider  a  family  of  dielectrics  for  which  the  electromagnetic  equations  can  be  solved.  The  param¬ 
eters  of  this  family  are  then  adjusted  until  the  scattering  coefficients  of  the  exactly  soluble  prob¬ 
lem  are  as  close  as  possible  to  the  prescribed  scattering  coefficients.  This  procedure  is  always 
approximate.  The  degree  of  success,  moreover,  depends  upon  one's  cleverness  in  choosing  the 
family  of  dielectrics. 

By  contrast,  the  procedure  for  finding  the  dielectrics  in  the  present  report  is  exact,  at  least 
in  principle.  If  computing  machines  are  available,  one  can  approximate  the  dielectrics  as  closeh 
as  one  wishes  even  when  one  cannot  solve  tlie  Gel'fand-l.fevitan  equation  exactly. 

This  report  contains  no  very  practical  solutions  because  the  technique  which  we  propose  is 
very  new  (it  was  developed  within  the  pa.st  ten  years),  and  most  of  the  work  in  the  past  has  been 
concern  ed  with  more  mathematical  aspects  of  the  problem.  It  is  hoped  that  this  report  will  stim¬ 
ulate  research  for  electromagnetic  applications.  Indeed,  its  sole  purpose  is  lo  stimulate,  if 
possible,  research  which  leads  to  practical  applications  of  the  theory. 
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II.  SCATTERING  AND  THE  ONE-DIMENSIONAL  SCHRODINGER  EQUATION 
The  time-dependent  one -dimensional  Schrddinger  equation  is 

^2 


T  k 


3 

3x^ 


+  V(x)  ^0(x,t) 


We  require  that  the  solution  i;i(x,t)  be  in  Hilbert  space  and  be  quadraticaily  integrable: 

r  k(x,t)i^ 

_oo 


dx  <  “ 


(1) 


(2) 


First,  let  us  assume  that  the  potential  function  V(x)  is  bounded  everywhere  and  V(x)  -►  0  as 
|x|  .  Later,  we  shall  make  V{x)  =  0  for  |x|  >  a.  We  proceed  to  solve  the  Schrddinger  equa¬ 

tion  in  the  usual  manner  by  separation  of  variables. 

Let  the  general  solution  be 


^00 

\  X(x|p) 


2  -iE.t 

g(p)  e'^-P  dp  +  X(x|E.)  g.  e  ^ 


(3) 


and  choose  x(x|p)  and  x(xiE^)  so  that  they  satisfy  the  follov/ing  equations,  respectively. 


— ^  +  V(x)  X(x|p)  =  X(x|p)  for  p^  >  0 
dx 


(4) 


+V(x)|  x(x|E.)  =  E.  x(x|E.)  forE.  <0 
dx 


(5) 


In  general,  the  solutions  of  these  second-order  differential  equations  are  not  unique  unless  bound¬ 
ary  conditions  are  imposed. 

The  discrete  values  E.  are  called  "point  eigenvalues,"  and  the  theory  of  Hilbert  space  re¬ 
quires  that  they  be  chosen  so  that 


y  X*(xlE.)  x(x|E^)  =  A.6.^ 


(6) 


where  A^  >  0. 

With  this  "boundary  condition,"  the  x(x|Ej)  become  unique  except  for  a  normalization  con¬ 
stant.  When  the  point  eigenvalues  E.  are  chosen  so  that  (6)  is  satisfied,  the  corresponding  x(x|Ej^) 
are  called  "proper  eigenfunctions."  Next  we  need  boundary  conditions  which  will  determine  the 
X(x|p)  uniquely.  Since  we  are  interested  in  solving  scattering  problems,  we  wish  to  impose 
boundary  conditions  such  that  the  resulting  set  of  x(x|p)  will  give  a  "simple"  description  of  scat¬ 
tering  problems. 

In  what  follows,  1st  us  assume  that  V(x)  is  bounded  everywhere  and  that  V(x)  =  0  for  |x|  >  a. 
Let  us  take  the  following  boundary  conditions: 
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For  p  >  0, 


X(x|p)  =  ^  ^  for  X  <  -a 

\l  Z-n  4^ 


J((x|p)  = for  X  >  a 
'^Ttt 


For  p  <  0, 

X(x|p)  =  ^  for  X  >  a 

42^  '/2t 


v(x|p)  = -^  forx<-a 

4^ 

Let  us  solve  the  specific  initial  value  problem  for  which 
=  0  for  all  i  , 


(7) 


(8) 


g(p)  =  0  for  all  p  <  0 


Then  from  (7)  and  (3)  we  obtain  a  particular  solution 
ipx  ._2 


yx.ii  = 


g,p,  dp  t  f 

o  4^  ' 


®  b(p)  g(p)  e"^^^  dp  forx<-a 


o  4^ 


0p(x,t)  =  j 


r“  eiP’^ _ -ip"t 


o  42ii 


t(p)  g{p)  e  dp 


for  X  >  a 


(9) 


For  any  fixed  p,  the  first  integral  in  the  expression  for  for  x  <  -a  and  the  integral  in  the 

expression -for  i/'  (x,t)  for  x  >  a  reduce  to  expressions  representing  waves  moving  in  the  positive 
x-direction.  The  second  integral  in  the  expression  for  i|ip(x,t)  for  x  <  -a  reduces  to  an  expression 
representing  a  wave  moving  in  the  negative  x-direction.  We  may  think  of  the  integrals  as  super¬ 
positions  of  such  waves.  With  this  interpretation  in  mind,  the  b(p)  are  called  the  reflection  coef¬ 
ficients  for  p  >  0,  and  the  t(p)  are  the  transmission  coefficients  for  p  >  0.  Hence,  given  the 
boundary  conditions  (7)  we  now  have  defined  x{x|p)  uniquely  for  p  >  0. 

In  a  similar  manner,  we  can  solve  the  initial  value  problem  for  which 


g.  -  0  for  all  i 


g{p)  =  0  for  p  >  0 

and  we  define  x(x|p)  uniquely  for  p  <  0  with  the  help  of  the' boundary  conditions  (8).  The  r(p)  and 
s{p)  are  called  the  reflection  and  transmission  coefficients,  respectively,  for  p  <  0. 

Let  us  call  x{x|p)  defined  for  p  >  0,  x^.(x|p)  and  x(x|p)  defined  for  p  <  0,  x.(x|p).  By  analyt 
continuation  we  may  define  x+(x|p)  for  p  <  0  and  xjxjp)  for  p  >  0.  It  can  then  be  shown  that 


X+(x|-p)  =  X*(x|p)  . 

X_{x|-p)  =  xt{x|p)  X 

X+(x|-p)  X.(x|p) 
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Also,  by  analytic  continuation  it  can  be  shown  that 


and 


b(-p)  =  b*(p) 
t(-p)  =  t*(p) 
r(-p)  =  r*(p) 
s(-p)  =  s*(p) 


for  p  >  0 


for  p  <  0 


s(p)  =  t(p) 


for  all  p 


Two  other  relations  satisfied  by  the  reflection  and  transmission  coefficients  are 

|b(p)  +  lt(p)  1^  =  1  for  p  >  0 

|r(p)  I  ^  +  I  s(p)  I  ^  =  1  for  p  <  0 

These  seem  to  imply  conservation  of  energy.  Also, 
lim  b(p)  =  0  , 

I  P 

lim  t(p)  -  1 

I P 1"*“ 

in  the  upper  half  plane.  This  is  always  true.  When  V(x)  dies  down  very  rapidly,  it  can  be  shown 
in  addition  that  b(p)  has  poles  at  the  points  p  =  i.^— E^;  that  is,  b(p)  has  poles  corresponding  to 
every  point  eigenvalue  E^.  This,  however,  is  true  only  when  V(x)  dies  down  sufficiently  rapidly. 
Examples  to  illustrate  this  will  be  given  later. 

Once  the  x(x|E^)  and  ,t(x|p)  have  been  defined  with  the  help  of  the  boundary  conditions  (6), 

(7)  and  (8),  we  find  that  they  satisfy  the  completeness  and  orthonormality  conditions.  If  chosen 
differently,  the  x(x|E.)  and  x(x|p)  would  not  satisfy  these  conditions. 

The  orthonormality  conditions  are 


X^ixlp)  x(x|p')  dx  =  6(p-p’) 


poo 

j  X*(x|E.)  x(xlE.)  dx  =  A.6.. 
X’^‘(x|p)  x(x|E.)  dx  =  0 


(10) 


The  completeness  relationship  is 


r 


X^CxIp)  x(x'|p)  dp 


X*(x|E.)  x(x'1 
Li 


Ep 


A. 

1 


=  6(x  —  x') 


(11) 


The  satisfaction  of  these  conditions  leads  to  a  complete  analogue  with  Fourier  transform  theory. 

Suppose  a  function  f(x)  is  in  L^.  Then  as  a  result  of  the  completeness  theorem  we  may  ex¬ 
pand  t(x)  as 
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(12) 


=  \  X(x|p)  g(p)  tip  +  E  X(x|E.)  g. 

^00  i  J. 

where 

roO 

X*(x|p)  f(x)  dx  , 

-00 

pOJ 

A.g.  =  j  x’'{x|E.)  f(x)  dx  . 

I'leiiLe  wo  n.ay  \/i.  ite 

f(x)  <=i>  g(p);  g. 

Also  suppose  that 

f^(x)  <;=>  g^(p);  g.^ 

r^(x)  g^{p); 

Then 

/^OO  iO  oO 

\  r^“^(x)  f^(x)  dx  =  \  g^*(p)  g^(p)  dp  +  2  g/’''g-^A.  .  (13) 

-00  -CO 

I 

But  this  is  simply  Parseval's  theorem. 

In  the  solution  of  scattering  problems  we  are  interested  in  the  reflection  and  transmission 

2  2 

coefficients.  In  particular,  for  p  >  0,  |b(p)|  and  for  p  <  0,  |r(p)|  give  us  the  relative  prob- 

ability  that  a  particle  of  momentum  p  will  be  reflected.  Similarly,  for  p  >  0,  |t(p)|  and  for 
2 

p  <0,  |s(p)|  give  us  the  relative  probability  that  a  particle  of  momentum  p  will  be  transmitted. 
The  relations 

|b(p)|^  +  |t(p)|^  =  1  , 

|r(p)|^  +  ls(p)|^  =  1  , 

state  that  one  event  or  the  other  will  oeeiir.  To  determine  the  reflection  and  transmission  coef¬ 
ficients,  Eq.  (4)  must  be  solved.  This  ean  be  done  in  the  usual  manner  by  matching  boundary 
conditions  where  V(x)  vanishes. 

For  an  alternate  integral  equation  te.'hnique,  we  write 

X(x|p)=-^-y^  f  e'IPl  Ix-x'|  (i4) 

nTZtt  IP' 

This  expression  for  x(x|p)  satisfies  the  wave  equation  and  all  boundary  conditions. 


5 


For  example,  take  p  >  0.  Then  (14)  gives 


X(x|p) 


eipx  ie-ipx 


X(x|p)  =-^— 


4^ 
ipx  I 


poo  .  j 

\  e^P’^  V(x')  x(x'|p)  ilx' 


71^  £ 


-ipx'  V(x')  x(x' I p)  dx' 


The  expressions  (15)  may  be  written  in  the  form  (7)  if  we  choose 


b(p) 


"'"'if  pj^  V(x')  x(x'|p)  dx 


t(p)  =  1 


»-ipx' 


V(x')  x(x'lli)  dx 


By  analytic  continuation,  b(p)  and  t(p)  can  be  defined  foij  p  <  0. 
Similarly,  take  p  <  0.  Then  (14)  gives 


ipx 


„ipx' 


x(x1p)  = 


e^ 

^/T?r 


/!  I  £  e--'  V,. 


V(x')  x(x'| 


1  + 


X')  X(x 


The  expression  (17)  may  be  written  in  tl  e  form  (8) 

r(p) 


[5)  dx'  for  X  >  a 

I p)  dx'  for  X  <  -  a 
kf  we  choose 


s(p)  =  1  + 


-  .  /^OO 

1  1  \ 

2  P  J..O 


eipx'  v(jj')  x(x'1p)  dx 


e'ipx'  v(x')  x(x'lil)  dx' 


'  p  >  0, 


By  analytic  continuation,  r(p)  and  s(p)  can  be  defined  fo:: 

As  an  example  of  the  solution  of  the  direct  problem  of  the  one-dimensional  Schrbdinger  equa¬ 
tion,  let  us  consider  the  delta  function  potential  V(x)  =  lB6(x),  and  let  us  take  B  >  0  so  that  there 
are  no  point  eigenvalues.  Hence  g-  =  0  for  all  i.  From 


From  (18), 


From  (14), 


b(P>  =-</!  I  2Bx(01p) 

t(P>=  1-/0 

r(p)  =  I  2Bx(01p) 

s(p)  =  i  +  JI  ^  2Bx(0|p) 

X(0|p)  =  -  forp>0 

NZa  (p  +  iB) 


for  p  > 


for  p  "H  0 


for  X  <  —a 


for  X  >  a 


(15) 


for  p  >  0 


(16) 


(17) 


for  p  <  0 


(18) 


(16), 
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Hence, 


'n/Ztt  (p  -  iB) 


for  p  <  0 


b(p)  = 


t(p)  = 


r(p)  = 


s(p)  = 


-iB 


p  +  iB 
P 

p  +  iB 

iB  1 
p  —  iB 

P 

p  —  iB 


for  p  >  0 


for  p  <  0 


By  analytic  continuation, 


b(-p) 

t(-p) 

r’(-p) 

s(-p) 


-iB 

iB 

-  b=''  (p) 

-p  +  iB 

p  —  iB 

-P 

P 

=  t’'‘(p) 

-p  -I-  iB 

p  —  iB 

iB 

-iB 

=  r*(p) 

-p  —  iB 

p  1  iB 

-P 

P 

=  s’^(p) 

-p  -  iB 

p  +  iB 

for  p  >  0 


for  p  <  0 


Hence  s(p)  =  t(p)  for  all  p  and  r(p)  =  b(p)  for  all  p,  the  latter  being  so  only  because  of  the  sym¬ 
metry  of  the  delta  function  potential.  Also, 

2  2 

lb(p)l^  +  |t(p)l^  =-2^ - 2+^ - 2"^  forp>0  , 

p  +  B  p  +  B 

2  2 

lr(p)r  +  ls(p)|^  -  -2^^ - 2+ -2^—=^  forp<0  . 

p  1  B  p  +  B 

III.  EXAMPLES  OF  ELECTROMAGNETIC  PROBLEMS  WHICH  CAN  BE  MAPPED 
INTO  A  ONE -DIMENSIONAL  SCHRODINGER  EQUATION 


We  shall  next  give  three  examples  of  one -dimensional  electromagnetic  theory.  It  v/ill  be 
shown  that  with  an  appropriate  mapping  each  example  can  be  reduced  to  a  one-dimensional 
Schrodlnger  equation. 


A.  Reflection  and  Transmission  of  Light  of  a  Fixed  Frequency 
at  Varying  *ngles  of  Incidence  by  a  Dielectric  Slab 

Consider  a  medium  in  which  the  dielectric  constant  e  =  6(x)  is  a  function  of  one  variable  x 
only  and  is  independent  of  time.  In  order  to  make  the  discussion  concrete,  we  assume  that 
£(x)  =  1  for  X  <  a  and  x  >  b.  Many  of  the  results  will  still  hold  if  eix)  dies  down  sufficiently  rap¬ 
idly  outside  the  range  a  <  x  <  b. 

Consider  a  ray  of  light  impinging  on  this  material,  the  angle  of  incidence  being  a.  The  plane 
of  incidence  is  the  xy-plane;  the  z-axis  projects  out  of  the  paper  (Fig.  i).  We  are  interested  pri¬ 
marily  in  the  amount  of  energy  reflected  and  transmitted  by  the  material.  It  will  be  shown  that, 
for  a  ray  polarized  so  th.at  the  electric  field  is  parallel  to  the  z-axis,  the  electromagnetic  equation 
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Fig.  1.  Roy  picture  showing  reflection 
and  transmission  of  light  by  a  dielectric 
slob. 


can  be  mapped  into  tlic  Schroclingcr  equation  and  tlie  b(p)  and  t(p)  of  this  associated  Schi-ddinger 
equation  will  turn  out  to  be  essentially  the  same  as  the  electromagnetic  reflection  and  trans¬ 
mission  coefficients. 

Using  Gaussian  units,  Maxwell's  equations  in  a  source-free  region  are 

V  ■  U  0 

V  ■  B  ^  0  ,  . 

V  X  K  =  —  i 

—  —  c  9t  ' 

VXU.if  ,  (19) 

where  c  is  the  frce-space  velocity  of  light.  We  also  have  the  constitutive  equations 
B  =  pH  , 

U  =  €(x)  B  .  (20) 

',V  i-  ,  osumc  that  p  =  p^,  where  p^  is  the  permeability  of  free  space  and,  to  eliminate  time  depend¬ 
ence,  write 

B(x,t)  -  B(x) 
n(x,t)  n(x) 

The  bar  under  x  indicates  dependence  on  both  x  and  y,  while  x  without  a  bar,  such  as  in  e(x), 
indicates  dc])endencc  on  the  one  variable  x  alone.  Sub.stituting  in  (19)  fi'om  (20)  and  (21),  we 
obtain 


V  •  c(x)  E(x)  0 
V  •  ll(x)  =  0 


ip  O) 

V  X  E(x)  =  — ^  l](.x) 


—  ic (x)  oj  ,,,  , 
V  X  n(x)  - E(x) 


Proceeding  in  the  usual  manner, 


V  X  V  X  ji(x)  =  VV  ■  E(x)  -  V  E(x) 


P  e(x) 

V  X  H(x)  =  — — 2 -  E(x) 

c 
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But 


Hence, 


V  ■  e(x)  E(x)  =  E(x)  .  Vc(x)  +  e(x)  V  ■  E(x) 


V  .  E(x)  - 


and 


Z‘'E(x)  + 


li^c(x)  0) 


E(x)  =  -  'I 


E(x)  ■  Ve(x) 
c(x) 


Let  us  write 


E(x)  =  E  ^(x)  +  E|^(x) 

where 

E^(x)  -  (O.O.E.^'xj) 


is  the  component  of  the  electric  field  normal  to  the  plane  of  incidence,  and 


E^(x)  ^  (Fyx).  Ey(x).  0) 

is  the  component  of  the  electric  field  parallel  to  the  plane  of  incidence.  Ihen  satisfies  the 

equation 


V^E  (X)  +  n^(x)  k^E..(x)  =  0 
—  z  —  ^ 


(22) 


where 


n(x)  - 

The  boundary  conditions  we  wish  to  impose  on  E^(x)  are 


ik-x  „  ik'-x 

E,(x)^  t 


\  2ir 


,p  ik-x 

E^(x)  = 

^  4zir 


■41^ 


lor  X  ^  a 


for  x  '■  1) 


^here  k  and  k'  are  vectors  given  in  terms  of  their  components  by 


(23) 


k  =  (k  cos  a,  k  sina,0) 
k'  =  (-k  cos  a,  k  sina,0) 

Then  the  first  term  in  the  expression  for  E^(x)  for  x  >  a  represents  a  wave  moving  in  the  direction 
(cos  a,  sin  0,0).  We  may  consider  this  to  he  an  incident  wave.  The  second  term  represents  a 
wave  moving  in  the  direction  (-cosa,  sina,0),  i.e.,  in  the  direction  of  specular  reflection.  It 
represents  a  reflected  wave.  B  is  the  reflection  coefficient,  and  |b1  gives  the  relative  pro¬ 
portion  of  energy  reflected  at  the  frequency  u).  Similarly,  the  expression  for  E^(x)  for  x  >  b 
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represents  a  wave  moving  in  the  direction  (cos  a,  sina,0).  It  represents  a  wave  transmitted 
through  the  material.  T  is  the  transmission  eocffieient,  and  |t|^  gives  ihe  relative  proportion 
of  energy  transmitted  at  the  frequency  o). 

We  shall  now  show  that  with  the  boundary  eonditions  (23)  on  E^(x)  we  can  reduee  equation 
(22)  to  the  Schrodinger  equation  by  an  easy  mapping  and  identify  B  and  T  with  b(p)  and  t(p), 
respectively. 

Let  us  consider  a  fixed  frequeney  co  >  0.  To  separate  .spatial  variables  we  write 
TP  /  \  iksinofy 

E^(x)  =  u(x)  e  y  .  (25) 

Wo  also  define 


p  =  k  cos  a 


V(x)  -  — k^(n^  —  1)  =  — k^(c(x)  —  1] 


(26) 


Since  oj  is  fixed,  p  is  a  function  of  the  angle  ot.  Usually,  we  will  be  dealing  with  c(x)  ^  1,  and 
hence  V(x)  <  0. 

Substituting  in  Eq.  (22),  we  obtain 


dx 


+  V(x) 


u(x)  =  p  u(x) 


(27) 


where  p  >0,  but  tliis  is  just  the  time-independent  Sclirodingor  equation. 
Suppose  tlie  boundary  conditions  we  impose  on  u(x)  are 


u(x)  = 


^  ,  b(p) 


\ZTr 


for  X  <  a 


2n 


u(x)  =  MP) 

'JTn 


for  \  >  b 


(28) 


Then  from  (25)  and  (26)  wo  obtain 


gi(keosax+ksinu;y)  ^i(-kcos  ax+k  sinoy) 


M  2ir 


„ik- X  .  k' 

e  ^  b(p)  e- 


slZii 


•\1 2ir 


for  X  <  a 


E,(x) 


t(p)  pifk  cos  «x+k  Sinn  y) 
\2n 


,,  ,  ik- X 
=  ‘■(P)  e - 

4Zn 


for  X  >  b 


If  we  identify  B  with  b(p)  and  T  with  t(p),  these  are  just  the  boundary  conditions  we  wish  to  im¬ 
pose  on  E^()c).  Hence,  we  can  solve  this  particular  electromagnetic  problem  by  using  the  mapping 
(26)  to  find  the  associated  Schrodinger  equation  and  identifying  its  reflection  and  transmission 
coefficients  with  the  reflection  and  transmission  coefficients  of  the  electromagnetic  problem. 
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We  might  note  again  that  in  quantum  mechanics  p  represents  momentum  and  — <«  <  p  <  «> .  How¬ 
ever,  p  in  Eq.  (27)  as  defined  by  the  mapping  (26)  is  a  function  of  the  optical  angle  of  incidence 
a,  and  this  imposes  restrictions  on  the  range  of  p.  The  technique  is  valid  for  complex  q  and 
might  therefore  be  used  to  handle  more  complicated  wave  fronts. 


B.  A  Transmission  Line  Problem 


The  second  example  we  shall  give  is  a  transmission  line  problem.  It  will  differ  from  the 
above  example  which  dealt  with  a  fixed  frequency  and  varying  angles  of  incidence  in  that  now  we 
shall  consider  a  fixed  angle  for  varying  frequencies. 

Consider  a  transmission  line  witli  distributed  impedances  l,(z)  and  capacitances  C(z).  The 
transmission  line  equations  are 


aV(z,t) 

9z 

81(z,t) 


-L(z) 


8I(z.t)  1 


8t 


_  _  _n(^)  8  V(z,  t) 

9z  ■  '  '  8t 


(29) 


where  V(z,t)  and  l(z,t)  are  the  potential  and  tlie  current,  rcspcctiveh  .  For  the  special  case 
L=L^  , 

C  =  C^  , 


we  obtain  at  once 


8  V(z,t) 
8z2 


L  C 
o  o 


8  V(z,t) 
8t^ 


This  has  a  solution  of  the  form 

V(z,  t)  =  F(z  -  vt)  +  g(z  +  vt) 


where 


1 

V  = - 

JL  C 

V  o  o 

Apparently,  we  have  again  a  wave  propagation  problem  in  which  reflection  and  transmission  coef¬ 
ficients  are  of  importance.  We  propose  to  show  that  this  problem  also  can  be  reduced  to  the 
Schrddinger  equation. 

Let  us  write 


V(z,t)  =  V(z)  e 
l(z,t)  I(z)  e"^“‘  . 

As  above,  we  take  L  to  be  a  constant,  i.e.,  L  =  L^.  However,  for  a  <  z  <  b  we  take  C  to  be  a 
function  of  z.  For  z  <  a  and  z  >  b  we  take  C  to  be  a  constant  C  =  C^.  Then 

^  =  ic.LJ(z)  , 

=  ia)C(z)  V(z)  , 
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and  hence, 


+  a.'^L^C(z)  V(z)  =  0 
9z 


(30) 


2  2 

We  define  k  =  u)  and  impose  the  following  boundary  conditions  on  the  solution 


V(z) 


ikz 
e _ 


,  -ikz 
B(k)  e 

sTZn 


for  z  <  a 


V(z)  = 


T(k) 

\(Tir 


for  z  >  b 


(31) 


We  shall  sliow  tliat  if  we  can  solve  a  certain  Schrodinger  equa+ion  then  we  can  solve  (30)  subject 
to  the  boundary  conditions  (31).  For  convenience,  let  us  rewrite  (30)  in  the  form 

f  k^(z)  V(z)  =  0  ,  (32) 

az 


wliere 


Next  let 


e(z) 


C(z) 

C 


7)  (Z)  =  C  ^^"^(Z) 


and  introduce  a  new  independent  variable  x  related  to  z  by 


-2 


(z) 


(33) 


or 


>?^z) 


dx 

dz 


We  also  introduce  a  change  of  dependent  variable, 
U  -  7)V  . 


(34) 


(35) 


Then 


Also, 


dV  _  1  ^ _ IL 

dz  "  ri  dz  2  dz 

4 


^  ^11 _ 2  dT)  du 

71  ,  2  2  dz  dz 

'  dz  7) 


i  dS  ,  2^ 

*■  3  'dz' 


^  2 
7)  dz 


^  - 

2  du 

dz 

''  dx 

d  u 

4  d  u 

2 

■  4  2 

dz 

dx 

+  27)' 


du  ^ 
dx  dx 
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d?) 

dz 


2 

=  V 


dx 


d^rj  _  4  d^rj 

^  2  ■  '2 
dz  dx 


Substituting  in  (32),  we  obtain 


4  d^u  ^  ,  3 

V  — p  +  2?) 
dx 


du  d?; 
dx  dx 


(„2  42w.2  ^ 


dx 


)  (V 


dx' 


u 

2 

V 


d^7) 

•> 

dx*' 


+  2i]' 


^2 

'  dx' 


j  2u  ,  i 

+  -Q  (>) 


di  ,2  2  4  u 

)  +  k  j)  — 

dx  '  '  7J 


Fi'om  this  we  have 


^2  ,2  , 

d  u  u  d  7)  ,  ,  2  . 

— — - ^  f  k  u  =  0 

dx  '  ax 


We  write  this  as 


dx 


+  qlx) 


u(x)  =  k  u 


where 


q(x)  = 


i 

n 


(36) 


(37) 


Equation  (36)  is  the  Sohrodingor  equation  corresponding  to  Eq.  (30)  or  (32).  Note  that  q(x)  is  in¬ 
dependent  of  frequency.  Tncidentally,  the  mapping  is  possible  only  if  e(z)  is  continuous. 

We  shall  now  turn  our  attention  to  an  investigation  of  the  mapping  and  the  shape  that  q(x)  will 
have  for  a  realistically  chosen  c(z). 

A  first  observation  is  that  x  is  a  monotonic  function  of  z.  Also,  since  the  first-order  dif¬ 
ferential  equation  relating  x  to  z  allows  for  one  constant  of  integration,  we  may  put  x  =  A  when 
z  =  a;  then  z  =  b  will  correspond  to  some  value  of  x  —  say,  x  =  B.  The  mapping  is  not  unique 

since  A  can  be  chosen  arbitrarily,  but  a  change  in  A  merely  implies  a  shift  in  x. 

Consider  a  particular  £(z)  such  that  e(z)  is  some  function  of  z  for  a  <  z  <  b,  and  £(z)  =  1  for 
z  ^  a  and  z  >  b.  Then  we  want  to  determine  first  how  x  depends  on  z  for  z  <  a  and  z  >  b.  Since 
X  is  a  monotonie  function  of  z, 

z  <  a  corresponds  to  x  <  A, 
z  >  b  corresponds  to  x  >  B. 

From  (34),  we  have  at  once  that 

X  —  A  =  z  —  a  for  z  <  a,  i.e.,  x  <  A  , 

x  —  B  =  z  —  b  for  z  >  b,  i.e. ,  X  >  B 

Suppose  Fig.  2(a)  represents  a  typical  €(z).  It  is  essential  that  £(z)  be  continuous  for  the  mapping 
determined  by  (33),  (34)  and  (35)  to  be  possible.  However,  it  is  not  necessary  that  the  derivatives 
of  c(z)  be  continuous.  Given  an  £(z)  as  shown  in  Fig.  2(a),  7)(x)  would  have  a  shape  as  shown  in 
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|3-3T-B03I(.-<)| 


Fig.  2.  {(z)  as  a  function  af  z  and,  corresponding  to  this  efz), 
the  general  shapes  taken  by  r|(x);  dri(x)/dx  and  d^(x)/dx2  as 
functions  of  x. 


z  z 

Fig.  2(b).  Furthermore,  dtj/dx  and  d  t)/d.x  would  have  shapes  as  shown  in  Figs.  2(c)  and  (d), 

respectively.  Note  that  Figs.  2(b),  (c)  and  (d)  are  not  Die  exact  representations  of  t)(x),  drj/dx, 

2  2 

and  d  t)/dx  for  c(z)  as  shown  in  Fig.  2(a).  but  merely  are  intended  to  show  the  general  shape  of 

these  functions  for  the  given  c(z). 

2  2  2  2 
q(x)  =  (l/»)(x))  (d  T|/dx  )  will  have  a  shape  very  similar  to  tliat  of  d  rj/dx  ,  the  important 

tiling  lacing  tliat  it  will  have  two  maxima  and  one  minimum.  Hence,  if  we  wish  to  solve  our  trans¬ 
mission  line  problem  tlirough  its  associated  Schrodingcr  equation  for  a  realistic  e(z)  as  shown  in 
Fig.  2(a),  we  must  solve  the  Scltrodinger  equation  for  a  q(x)  such  as  shown  in  Fig.  3.  We  might 
note  that,  if  the  bumps  in  Fig.  3  were  replaced  by  delta  function  like  discontinuities  such  a  q(x) 
would  still  correspond  to  a  realistic  €(z).  Since  q(.x)  =  0  for  x  <  A  and  x  >  B,  we  again  have  a 
scattering  problem  and  shall  look  tor  solutions  of  the  Schrodingcr  equation  which  describe  scat¬ 
tering.  We  shall  sliow  tliat  if  wc  can  solve  this  scattering  problem,  we  will  also  have  solved  the 
transmission  line  problem  and  will  be  able  to  relate  the  coefficients  B(k)  and  T(k)  to  the  trans¬ 
mission  and  reflection  coefficients  of  the  Schrodingcr  equation.  Let  us  take  as  solutions  of  the 
Schrodingcr  equation 


u(.x) 

u(x) 


u/i  >  "iltx 
b(k)  e 

•^fzi 


...  .  ikx 
t(k)  e 


for  X  <  A 


for  X  >  B 


But  when  x  <  A  we  have  z  <  a;  17  =1  and  therefore  V(z)  =  u(x). 


Also,  .X  —  A  =  z  —  a. 


(38) 


Hci.ee 
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V(z)  = 


+  b(k)  gik(a-A) 

hjzi  ‘'Tzi 


gikz  g-ik(a-A) 


-ik(a-A) 


^  +  b(k)  32ik(a-A)l 

■  4^  'J  Zir  J 

_  g-ik(a-A)  ^  B(k) 

I-n/^  4^  J 


where 


2ik{a-A) 


B(k)  =  b{k)  e 

Similarly,  when  x  >  B  we  have  z  >  b;  7,  =1  and  therefore  V(z)  =  u(x).  Also,  now  x  -  B  =  z  -  b. 


Hence, 


Ni  Zlt 
^ikz 
4zir 


t(k)  e^^^  ^-ik{a-A)  ^-ik(b-B)+ik(a-A) 


_  T(k)  e^*^^  ^-ik{a-A) 

4ZTr 


where 


T{k)  =  t{k)  e 


ik{a-b+B-A) 


Fig.  3.  q(x)  as  a  function  of  x. 


The  factor  merely  changes  the  incident  wave.  Hence,  we  can  obtain  a  solution  of  our 

transmis.<-  nn  line  equation  (32)  subject  to  the  boundary  conditions  (31)  with  the  help  of  the  map¬ 
ping  given  by  (33),  (34)  and  (35).  We  solve  the  Schrodlnger  equation  (36)  subject  to  the  boundary 
conditions  (38)  and  obtain  the  B(k)  and  T(k)  of  the  transmission  line  problem  from  the  reflection 
coefficients  b(k)  and  the  transmission  coefficients. t(k)  of  :he  solution  of  the  Schrddinger  equation 

by  using  the  relationships 


B(k) 

T(k) 


b(k)  e 


2ik(a-A) 


t(k)  e 


ik(a-b+B-A) 


(39) 


In  practice,  reducing  the  transmission  line  equation  to  the  Schrddinger  equation  involves  the  fo^ 
lowing  steps.  Given  e(z),  we  use  (33)  and  (34)  to  find  x  as  a  function  of  z.  The  choice  of  the 
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constant  of  integration  A  is  arbitrary  but  is  usually  determined  by  symmetry  considerations. 
We  now  have  e  and  rj  as  functions  of  x  and  can  use  (37)  to  compute  q(x). 

Next,  we  shall  generalize  the  above  problem  slightly.  Suppose  we  again  take  the  transmis¬ 
sion  line  equation  (32)  and  let  e(z)  =  1  for  z  <  a;  c(z)  =  some  function  of  z  for  a  <  z  <  b;  but 
c(z)  =  Cp  =  constant  ^  1  for  z  >  b. 

We  shall  require  that  V(z)  satisfy  the  following  boundary  conditions; 


ikz 


V(z)  =■ 


B(k)  e 


ikz 


V(z)  = 


'•JZTt 


for  z  <  a 


for  z  >  b 


(40) 


The  mapping  is  independent  of  the  boundary  conditions,  and  we  shall  use  the  same  mapping  as 
above.  Again,  we  choose  z  =  a  to  correspond  to  x  =  A,  and  since  x  is  a  monotonic  function  of 
z  we  have  z  <  a  corresponds  to  x  <  A.  Also,  since  for  z  <  a  wo  choose  c(z)  -  1,  wo  have  as  be¬ 
fore  X  —  A  =  z  -  a  for  z  <  a  and  x  <  A.  z  =  b  will  correspond  to  some  value  of  x  —  say,  x  =  B  — 

1/ 

and  z  >  b  will  correspond  to  x  >  B.  However,  for  z  >  b  wo  now  have  c(z)  =  ^^d  hence  i;  = 

Therefore, 


and 


Again  lot  us  assume  a  simple  realistic  shape  for  e(z)  such  as  shown  in  Fig.  4(a).  The  apjjroxi- 

2  2 

mate  shapes  for  ij,  di;/dx  and  d  tj/dx  are  shown  in  Figs.  4(b),  (c)  and  (d).  q(,\)  will  of  course 

2  2 

again  have  a  shape  similar  to  that  of  d  p/dx  .  Therefore,  in  this  case  we  shall  be  dealing  with 
a  q(x)  having  a  shape  like  that  shown  in  Fig.  5, 

We  solve  the  Schrddinger  equation  for  this  q(x)  subject  to  the  boundary  conditions  (38).  For 
X  <  A  we  have  z  <  a;  i)  =  1;  V(z)  =  u(x);  and  x  -  A  =  z  -  a  as  before.  The  boundary  conditions  on 
V(z)  tor  z  <  a  are  also  the  same  as  above.  Hence,  again  we  obtain 


B(k)  --  b(k)  . 

For  X  >  B  we  have  z  >  b;  i)  =  V(z)  - 


and  X  —  B 


Hence, 


V(z) 


t(k)  iMB-TT^b) 


-lk(a-A|  IW-Tr^btB-AI 

In  this  case,  the  transmission  coefficient  t(k)  of  the  associated  Schrbdinger  equation  and  the  T(k) 
of  the  transmission  line  equation  are  related  in  the  following  manner; 


17 


T(k) 


t(k) 


1/4 


and  again  we  are  able  to  solve  the  transmission  line  problem  by  mapping  it  into  the  Schrbdinger 
equation. 

C.  Reflection  and  Transmission,  for  Varying  Frequencies, 
of  a  Plane  Wave  Normally  Incident  on  a  Dielectric  Slab 

Finally,  we  shall  show  that  Eq.  (32),  the  transmission  line  equation,  also  occurs  in  a  simple 
electromagnetic  problem.  In  this  case,  k  and  e  are  defined  differently,  of  course. 


I  i-3T-8035 


Fig.  6.  Plane  wave  incident  normally  an  a  slab 
of  dielectric  constant,  e(z). 

Consider  normal  incidence  of  a  plane  wave  traveling  in  the  z-direction  on  a  slab  of  material 
whose  dielectric  constant  e  =  e(z)  varies  in  the  z-direction  (Fig.'b).  It  was  shown  previously 
that,  assuming  a  time  dependence  e'^‘^*'  and  using  Ma-xwell's  equations  and  the  constitutive  equa¬ 
tions,  we  find  that  the  electric  field  E(x)  must  satisfy  the  equation 

2 


+  -2—^ - 


E(x)  =  -  V 


rE(x)  •  ve(z) . 
I  TIF)  ■ 


(41) 


Since  we  are  now  assuming  that  e  =  €(z)  varies  in  the  direction  of  propagation  and  since  E(x)  must 
be  normal  to  this  direction,  we  see  at  once  that  the  right-hand  side  of  (41)  vanishes. 

Let  us  assume  further  that 

E(x)  =  [E^{z),0,0]  . 

Then  we  have  at  once 


9^E  (z)  p  u>^ 
X  ,  o 
_ 


^  e(z)  E^(z)  =  0 


9z  c 

Similarly,  if  we  had  assumed  that 

E(x)  =  [0,Ey(z),0]  , 

we  would  have  obtained  the  same  equation  for  Ey(z),  namely. 


9^E^ 


9z 


^  e(z)  Ey(z)  =  0 


18 


We  can  write  cither  of  these  in  the  form  of  the  transmission  line  equation  (32): 


3  E  (z) 


8z 

where  now  v  =  x  or  y. 


2  +  k  e(z)  E^(z)  =  0 


(42) 


u)  [Ji  e 

,2  '^0  0 


e(z)  = 


c 

f(z) 


and  e  is  the  dielectric  constant  of  free  space., 
o 

IV  THE  INVERSE  PROBLEM  FOR  THE  SCHR5dINGER  EQUATION  - 
THE  GEL'FAND-LEVITAN  ALGORITHM  AND  ITS  PROOF 

The  direct  problem  of  solving  the  two  time-independent  equations  arising  from  the  Schrodinger 
equation  has  been  treated  earlier  in  this  report. 

To  recapitulate  briefly  in  this  case  specifying 

(a)  V(x), 

(b)  A^  the  normalization  of  x(-xlE^): 

A.  =  f  lx(x|Edl2dx  . 

^  -00 

(c)  The  boundary  conditions  on  x(xlp); 

lim  x(x1p)  = 


+  b(p)  e'^P^ 


V2¥ 


lin.  X(xlp)=lii^ 

v-«o  WZTT 


for  p  >  0 


e^P^  ,  r(p)  e-'P^ 

lim  x(x  I  p)  =  - 

x-«o  w  Ztt  \f  Z  tt 


lim  X(xlp)=-^M:|J^ 
4ZT! 


for  p  <  0 


we  solve  the  two  equations 

,2 


d 

dx^ 


+  V(x)|  x(x1p)  =  P^X(xlp)  for  p  >  0  , 


+  V(x)]  x(x I  E.)  =  E.x(x  1  E.)  for  E .  <  0  , 

dx 


and  determine  E.,  b(p).  t(p),  r(p).  s(p),  x(x  I  P)  and  x(x  1  E.) .  It  was  also  shown  that  the  x(xl  p)  and 
X(xlE.)  thus  determined  satisfy  the  completeness  and  orthonormality  conditions,  and  hence,  that 
any  function  f(x)  can  be  expanded  in  the  following  manner : 
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where 


f(x)  =  I  x(x|p)  g(p)  dp  +  2  X(x|E.)  g.  , 


g(p)  = 

A.g.  = 
1^1 


X*(x|p)  f(x)  dx 


X*(x|E.)  f(x)  dx 


We  shall  next  discuss  the  inverse  problem,  namely,  the  problem  of  determining  V(x)  from  a 
knowledge  of  A.,  E.,  b(p)  and  the  boundary  conditions  on  x(x|p).  This  problem  was  first  treated 
by  Gel'fand  and  Levitan.  However,  they  dealt  with  the  radial  equation  for  f  =  0,  and  we  shall  be 
dealing  with  the  one -dimensional  equation  for  the  range  -»  <  x  <  «> . 

We  shall  first  give  the  Gel'fand -Levitan  algorithm  and  then  proceed  to  the  proof.  The  proof 
will  be  given  in  two  parts. 

Gel'fand-Levitan  Algorithm 
If,  given  b(p),  Ej^,  and  A^,  we  define 


R(x) 


^  b(p)  e“^P^  dp  +  ^ 


-E.x 
1 


A; 


(43) 


and  assume  that  the  Gel'fand-Levitan  equation  for  x  >  y, 
K(x,y)  =  -R(x  +  y)  -  r  K(x, z)  R(z  +  y)  dz 

—Ob 


(44) 


has  a  unique  solution  for  K(x,y),  then  the  potential  V(x)  is  given  by 
V(x)  =  2^  [K(x,x)]  , 

and  for  p  >  0 , 

^  +  -j-  r  K(x,y)  [eW  4  b(p)  e-W)  dy 

’sTZt!  “/zir 


/-E.x  nx 

X(x|E.)  =  e  +  \  K(x,y)  e  dy 


(45) 


(46) 


In  principle,  the  inverse  problem  which  requires  the  solution  of  a  linear  integral  equation  — 

the  Gel'fand-Levitan  equation  —  is  no  more  difficult  than  the  direct  problem  which  requires  the 

solution  of  a  linear  differential  equation.  We  note  that  for  fixed  x  the  Gel'fand-Levitan  equation 

is  a  Fredholm  equation  and  hence  is  of  a  well  understood  standard  form, 

Also,  as  x--"  we  have  R(x)  0.  Each  of  the  terms  of  the  sum  S  e^  i  /A^  certainly  tends 

i 

to  zero  as  x”**  and  with  the  help  of  the  Riemann-Lebesque  lemma  it  can  be  shown  that  the 
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integral  in  the  expression  for  R(x)  vanishes  as  x  -  -«>.t  The  given  b(p)  will  of  course  satisfy  the 
relations  b(-p)  =  b*(p),  lim  b(p)  =  0  in  the  upper  half  plane.  We  might  point  out  that  as  a 

IpK°° 

consequence  of  the  first  of  these  two  relations  R(x)  is  real  and  this  in  turn  implies  that  K(x,y) 
is  real. 

In  general,  it  is  necessary  that  b(p),  E^,  and  are  specified  before  the  Gel'fand-Levitan 
algorithm  can  be  applied.  However,  for  very  short  range  potentials  b(p)  and  E.  are  sufficient 
and  need  not  be  given. 

Let  us  specify  two  functions  i/)(x|p)  and  i/)(x|e^)  which  are  solutions  of  the  equations 


+  V(x)]  !/)(x1p)  =  p^!/<(x1p)  for  p^  >  0 


dx 

,2 


--^+V(x)  i/)(x|E.)  =  E.i/)(x|E.)  forE.  <0 

dx 


(47) 


The  !/<(x1e.)  are  to  be  identical  with  the  x(x|E.)  and  hence  are  determined  by  imposing  the 
same  boundary  condition  on  them  as  was  imposed  on  the  x(x1e.),  namely, 


y  i/<*(xlE.)  i/-(xlE.) 


dx  =  A.a.. 


for  all  p 


(48) 


However,  the  boundary  condition 

,  c^P"" 

lim  !/'(x|p)  =  —ir 

which  we  impose  on  the  !/»(x|p)  is  different  from  tlie  condition  imposed  on  x(x|p).  One  reason  for 
the  choice  of  this  boundary  condition  is  convenience.  Also,  with  this  condition  the  !(i(xlp)  are 
analytic  functions  for  all  p.  The  two  independent  solutions  i(i(x|p)  and  !/)(xl-p)  may  now  be  super¬ 
imposed  and  we  write 

X(x|p)  =  !(i(x|p)  +b(p)  !/'(x|-p)  for  p  >  0| 

X(x|p)  =  !/)(x1p)  +  r(p)  !/>(xl-p)  forp<0 
Next  we  proceed  to  Part  I  of  the  proof  of  the  Gel'fand-Levitan  algorithm.  In  this  part  of  the  proof 
we  shall  show  that 


(49) 


^(x1p)=—  +  —  r  K(x,y)e'Pydy 


i/)(x|  E 


;.)  =  '  +  \  K(x,y)  e^  dy 

1 


(50) 


satisfy  the  equations  (47)  when  (45)  is  satisfied.  It  is  immediately  apparent  that  !f(xlp)  as  given 
by  (50)  satisfies  the  boundary  condition(48). 


tThe  Riemonn-Lebesque  theorem  (Lighthill,  Introduction  to  Fourier  Angiysis  ond  Gerierolized  Funetj^s,  p.^) 
stotes  thot:  If  b(p)  is  on  ordinary  function  obsalutely  integroble  from-oo  to  co,  then  its  Fourier  transform  -  0  os 
I  x|  -  ®.  However,  b(p)  in  this  cose  consists  of  a  given  set  of  values  ond  is  therefore  certomly  obsolutely  m- 

tegroble  over  {— <»,oo). 
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Next  we  calculate  (d^/dx^)  ;/'(x|p)  +p^;/)(x|p). 


dx 


^(xip)  =  iP£^  +  'SM  r 

'  '  u. 


jipx 

'v/TtF  ’sfZn 


\IZti 

2  e^P^  ipK(x,x)  ipx  ,  e^P^  d 


iyx,y)  e'Py  dy 


dx 


^,(x|p)  =  _p-  2:^  +  e--  ^  K(x,x) 

^  ^/2^  4Zn  Vli 


K  (x,x) 

+  e^P^  +  ■ 


=r 

^TT 


Kxx(>^-y)  e^'^^dy  , 


where 


iyx,x)  =  ^  [K(x,y)]|y^^ 


Ky(x,x)  =  ^  [K(x,y)]|y^^ 


(51) 


Integrating  by  parts  and  making  use  of  the  fact  that  as  x  R(x)  -*  0,  and  therefore  also 

K(x,y)  —  0,  we  have 


,  2  „ipx 

p^^{x\p)  =  £ -  + 


'7^  _ 


n/TtT  ^  -00 

p^  e^P^  _  ipK(x, x)  ipx 


K(^,y)  p^  e^P^  dy 


=r 

Zir  *^-=0 


e  '  "  + 

\!^  \IZTr  ’sTUr 


2  ipx  „v  ; _  I\,(X,X) 


K  (x,y)  ip  e^P^  dy 


_  p"  e  _  ipK(x,x)  gipx  ^  y 
^fzii  's/T?  "^Zir 


^ipx 


-^r  K 

^/2i  yy 


(x,y)  e^P^  dy 


Hence, 


+p2^,(x|p)  =-^  [K(x,x))  +K  (x,x)  +K  (x,x)} 


dx 


'4Zv 


C  [K  (x,y)-K  (x,y)]e^py 


^fZir 

7  p^P^  d  1 

4  [K(x,x)]  +—  [K  ' 
V27  -slzi 


dy 


[K^v(x,y)  -  Kyy(^-y)l  dy  .  (52) 


Since 


K^(x,x)  +  KJx,x)  =  [K(x,x)]  , 

C  ( 
yy 

K(x,y)  =  -R(x  +  y)  -  C  K(x,  z)  R(z  +  y)  dz  . 

v  _QQ 


-X -  -  y'  '  ■  dx 

we  proceed  to  find  an  expression  for  ~  have 


(53) 
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Henee, 


K  (x,y)  =  -R'(x  +  y)  -  K(x,x)  R(x  +  y)  -  C  K  (x,  z)  R(z  +  y)  dz 
K^j^(x,y)  =  -R"(x  +  y)  -  K(x,x)  R'(x  +  y)  -  ^  [K(x,x)]  R(x  +  y) 
-K^(x,x)  R(x  +y)  -  j'  K^^(x,z)  R(z  +  y)  dz  , 


K  ,(x,y)  =  -R'(x  +  y) 


-r 


z)  R'(z  +  y)  dz 


=  -R'(x  i  y)  -  K(x,  x)  R(x  +  y)  f  [  K  (x,  z)  R(z  +  y)  dz 

^  -OO 

K  (x,y)  =  -R"(x  +  y)  -  K(x,x)  R'(x  +  y)  +  y  K^(x,  z)  R'(z  +  y)  dz 


-R"(x  +y)-K(x.x)  R'(x  f  y)  +K  (x,x)  R(x  +  y)  -  C 

y  ^  _« 


K^^(x,  z)  R(z  +  y)  dz 


K^^(x,y)  -  Kyy(x,y)  = [K(x,x)]  +K^(x,x)  +Ky(x,x)}  R(x  +  y) 


Of  using  (53), 


’^xx'^^'y’  "  ^  (l<(x,x)]  R(x  +  y) 


This  is  an  integral  equation  for  Ix^^^lx.y)  —  R^^,(x,y).  Let  us  try  to  find  a  solution  in  the  forr 

Kxx(x.y)  -  Kyy(x,y)  =  L(x,y)  2  £  [K(x,x)]  .  I 

Substituting  (65)  in  (54),  wo  obtain 


L(x,y)  ^  -R(x  +  y)  -  C  L(x,  z)  R(z  ^  y)  dz 

''  _oo 


Hence,  L(x,y)  satisfies  the  same  equation  as  K(x,y).  If  we  assume  that  this  equation  has  a  unique 
solution,  then  L(x,y)  =  K(x,y)  and 

Kx^(x,y)-Kyy(x,y)  =  2K(x,y)  ^  K[x,x]  .  (57) 
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Using  (57)  in  expression  (52),  we  obtain 


dVjp)  +p2^(,|p)  =  i^^[K(x,x)] 
dx 

+  —  r  2K(x,y)  ^  [K(x,x)]  e'P^  dy 
sfZir 

^  2d[K^.^  [s!!!  +  _L  f  K(x,y)  e^py  dyl 

^  42^  J-oo  I 

=  2  [K(x,x)]  0(x|p) 

Hence,  if 

V(x)  =  2  ^  [K(x,x)]  , 

then  i;’(x|p)  satisfies  Eq.(47).  The  proof  for  0(x|Ep  is  identical.  This  completes  Part  1  of  the 
proof. 

In  Part  1  of  the  proof  of  the  Gel 'fand -Levitan  algorithm  we  have  shown  that  the  i/>(x  |  p)  and 
i)i(x|E.)  as  given  by  (50)  satisfy  the  equations  (47).  Also,  ^Ixlp)  obviously  obeys  the  boundary 
condition  (48).  In  Part  II  of  the  proof  we  shall  prove  the  following  completeness  relationship; 


;;i(x|p)  >p*(x' 


;/i(x|p)  b(-p)  ;;i*(x'|-p)  dp 


i 


;;i(x|  E.)  i/'(x'|  E. 


(5(x  -  x') 


(58) 


We  shall  then  use  this  completeness  relationship  to  show  that  the  t''(x|p)  arc  linearly  independent, 
that  the  ;/)(x|E.)  are  quadratically  integrable,  have  the  normalization  A^  and  are  therefore  the 
point  eigenvalues  of  the  Schrodinger  equation  (47),  and  finally,  that  the  b(p)  as  used  in  (43)  are 
indeed  the  reflection  coefficients  of  the  fichrddinger  equation. 


Let  us  first  prove  (58).  For  x  >  y  we  have 

K(x,y)  = -R(x  +  y)  -  C  K(x,  z)  R(z  +  y)  dz 

aO 


=  -R(x  +  y)  -  Y 


K(x,  z)  b(p)  e  dp 


J^Az+y) 

_r 

>^-oo  .  1 


=  -R(x  +  y) - ^ 

42An 


0(x|-p) 


^I^]b(p)  e'^Py  dp 
4Zn 
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Therefore, 


K(x,y)  = - -  \  ^(x|-p)  b(p)  dp  -  V  - — 

kI  7‘rr  ^  _Qrt  ^  i 


A. 

1 


(59) 


Also, 


K(x,y)=-r  ^l>{x|-p)  b(p)  f^(>(y|-p) - L.  K(y,z)  e’^P^  dz 


dp 


E 


!/<(x|E.) 

~A. 

1 


ny\E.)-^' 


y 

K(y,z)e^  ^  dz 


/^OO 

=  -\  i;i(xl-p)  b(p)  ^)-(y|-p)  dp  -  ^ 

.-o<i 


Wx|  E.)  i/>(y  I  E.) 


A. 


f*y  1 

+  \  K(y,2)  — \  ^)i(x  |-p)  b(p)  e"^P^  dp  +  y 

^ -oo  ^fZir  -oo 


/-E.zt 
if(x|E.)  e*^  ^  ' 


A. 


dz 


But,  according  to  (59),  the  expression  in  brackets  in  the  last  integral  is  just  -K(x,z).  Hence, 
rearranging  terms  v;e  now  have 


/-»0O 

\  i/'ixl-p)  b(p)  i)<(y|-p)  dp  t  Y, 

_QO 


,  (/-(x|E.)  Hy\E.) 

— X - ^  =  -K.x,y) 


-r 


K(y,z)  K(x,z)  dz 


(60) 


This  is  true  for  x  >  y.  However,  the  left-hand  side  of  (60)  is  symmetric  in  x  and  y  so  that  for 
X  <  y  we  liave 


\  i/'(y|-p)  b(p)  i/'(y|-p)  dp  +  ^ 

-Ort 


«y|E.)  ^/.(x|E,) 


A. 

1 


■K(y,x) 


-r 


K(x,  z)  K(y .  z)  dz 


Introducing  the  Heaviside  step  function, 

1  when  X  >  0 

llx)  = 

0  when  x  <  0  ,  (61) 

and  replacing  p  by  -p  in  the  integral  on  the  left-hand  side,  we  finally  have 

#(x|E  )  ^/.(y|E  ) 

\  4>{x\p)  b(-p)  i/'  '‘(y|-p)  dp  +  2,  - ^  "  -T](x-y)  K(x,y)-  p(y  -  x)  K(y,x) 

-00  .  i 

1 

-9(x-y)C^  K(y,z)  K(x,z)  dz -T)(y  -  x)  C  K(x,  z)  K(y,  z)  dz  .  (62) 

_oo  -OO 
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Next,  using  (50),  we  write 

pcO  •  noO 

\  '^’{xIp)  fUy\p)  dp  =  \ 

^-00  '-'-oc 


^ipx 


Ztt  'slZir 


=r 


K(x,  z)  e*^"  dz 


But 


/Vlso , 


.  1  P  . . „-ip» 


“V 

nz  J_. 


K(y, v)  e  ^  dv 


dp 


sfZTT  ^FZtI 

=  ^  dp  +  2^  K(y,v)  dv  d, 

t  2*^  y  K(x,z)  dz  ^  dp 

+  ^  r  KK-,z)dz  r  K(y,v)dv  f  dp  . 

-00  '^-OO 


i‘°0 

gip(x-y)  Zirdix  -  y) 


^  r 

27r  J 

^  .or 


K(y,v)  dv  C  dp  =  K(y,v)  a(x  -  v)  dv 

^  -uO 


=  y  K(y,> 

-m3 


Ztt 


f  K,x,z)dzy 


,ip('^-y) 


K(y,x)  whun  X  <  y 

0  when  x  >  y 

r)(y  -  x)  K(y,x) 

pX 

dp  \  K(x,  z)  6(z  -  y)  dz 

— oo 


K(x,y)  when  y  <  x 

0  when  y  >  x 

))(x  -  y)  K(x,y) 


r  K(x,z)dzy  K(y,v)dvr  e'P‘''^-^>  dp  -  r  K(x,z)dzr  K(y,v)dva(v-z) 

J_0O  -^-00  -^-00 


Ztt 


f'y 

^  \  K(x,z)  K(y,z)  dz  when  x  >  y 

-CO 

=  r  K(x,  z)  K(y,  z)  dz  when  y  >  x 

^  -oO 

=  ;,(x-y)y  K(x,  z)  K(y,  z)  dz  +  T)(y  -  X)  ^  K(x,  z)  K(y ,  z)  dz 
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Hence, 


^oO 

\  !;'(x1p)  i/»’’‘(y  Ip)  dp  =  6(x  -  y)  +  riiy  -  x)  K(y,x)  +  t)(x  -  y)  K(x,y) 

v  -Oft 


+  7!(x 


K(x,  z)  K(y,  z)  dz  +  tj(y  —  x)  T  K(y,  z)  K(x,  z)  dz 

d  _ocj 

Adding  (6E)  and  (63),  we  obtain 

poo  poo 

\  i/'(x|p)  0’'(ylp)  dp  +  \  ifCxlp)  b(-p)  (('’’“(yl-p)  dp 

^  -QO  ^ -to 

!))(x|E.)  i/'(ylE.) 


z 


A. 


--  6(x  -  y) 


(63) 


Tills  completes  the  proof  of  the  completeness  relationship  (58). 

A  direct  consequence  of  (58)  is  tliat  any  quadratically  integrable  function  t(x)  can  be  expanded 

in  the  following  manner: 


f(x)  ^  I  ^'(xlp)  a(p)  dp  +  Y.  a. 

-QO 


(64) 


where 


all.)  =  £ 


f(y)  I/'* (yip)  dy  +  \ 

^  -OC 

f(y)  <,''(ylE.) 


((y)  b(-p)  tl'’“(yl-p)  dy 


A, 


dy 


(65) 

(66) 


Let  us  treat  two  special  cases; 

(a)  t(x)  =  t''(x|p) 

(b)  f(x)  =  0(x|E.)  . 

0(x1p)  is  not  quadratically  integrable,  but  it  is  still  symbolically  possible  to  treat  case  (a), 
(a)  When  f(x)  -  0(xlp),  we  see  at  once  that  in  this  case  we  must  have 
a(p)  =  6(p  -  p') 


a.  =  0 
1 


The  first  of  these  gives 


i/<(ylp')  I'^'lylp)  dy  +  ^  ./.(y  1  p')  b(-p)  (yl-p)  dy  =  6(p  -  P')  •  (67) 


The  second  gives 


£ 


iflylp')  '/'(yl E.) 


A. 

1 


-  dy  =  0 


(68) 


Hence,  the  i/»(xlp)  and  0(xlEi)  are  orthogonal. 
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(b)  When  f(x)  =  ^(xIe.),  we  see  at  once  that  in  this  case  we  must  have 


a(p)  0 

a.  =  6.. 

1  1] 

From  the  second  of  tlicse  we  obtain 
ii)(ylE  )  i/>(y|Ep 


r 


A. 


dy  =  6.. 
^  ij 


^  vMyjE.)  ,(y|K.)  dy  -  A.6.. 


Hence,  the  i;i(x]E^)  are  qiiadralically  integrable  and  have  the  proper  noi-malization. 

We  shall  now  prove  that  the  b(p)  as  used  in  the  Gcl'fand-Levitan  algorithm  are  indeed  tlie 
reflection  coefficients . 

l,ct  us  consider  the  outgoing  wave,  i.e.,  \(x|p)  for  p  >  0,  and  assume  that  the  reflection  coef¬ 
ficient  is  b(p)  and  not  b(p).  With  the  help  of  the  completeness  relationships  for  tf<(x|p)  and  x(\lp), 

/\ 

we  sliall  sliow  tliat  b(p)  e  b(p). 

Since  x(>tjl))  and  I|))  are  solutions  of  the  same  Schrodingcr  equation  foi'  difforciit  boundary 
conditions,  we  may  write  x(-'<l|>)  as  a  linear  combination  of  the  'wo  i’Hlependcn'  solutions  i!'(s|p) 
and  f/<(xl-p)-  Wo  write 

X(x|p)  i|(x|p)  +  b(p)  tMxl-p)  for  p  >  0  (70) 

We  liavc  cliosen  i('(x|p)  so  tliat  it  satisfies  tlie  boundary  condition  (-IH),  and  as  a  result  of  this  x'l^^lp) 
as  given  by  (70)  will  satisfy  the  correct  boundary  condition 


lim  x(>;lp) 


\.T? 


h(p)£‘'^^ 

\  2- 


for  p  >  0 


However,  the  x(>'iit)  satisfx'ing  this  boundary  condition  is  uniciue,  and  hcnct?  b(|))  is  the  reflection 
coefficient  as  defined  previously. 

Similarly,  unu|ueness  of  a  solution  for  a  given  boundary  condition  gives  us 

V(x|  p)  =  t(p)  i|(xl|))  for  p  <  0  .  (71) 


A  ^  2.  ^  I- 

t(p)  is  the  transmission  coefficient ,  and  |b(p)l  +  lt(p)|  -  '  lor  all  p. 

Again  we  clioose  the  proper  eigenfunctions  x(x|E^)  and  Clx]  ICp  to  be  identical  and  to  have  tlie 

same  normalization  A^.  We  might  point  out  that  since  K(x,y)  i.:.  real,  ij'(x|  b^)  and  x(-x|l''j)  are  real 

We  now  liave  tlie  completeness  relationship  (H)  for  the  x(xlp)  and  x(x|  namely, 

X(x]E.)  x(-'^'|1^^) 

\  X(x|p)  X*(x'|p)  'Ip  i-  L  - A -  =6(x-x') 

d_0O  .  i 

1 


and  the  completeness 


relationship  (58)  for  the  ^i(x|p)  and  i/i(x|Ej^)  =  x(xlr'q),  namely, 

poo  ‘(’(xlE.)  i/dx'lE.) 

Ip)  up  +  \  tf'(x|p)  b(-p)  i/i^fx'l-p)  dp  +  2.,  - ^  = 

— OO  •  1 


6(x  —  x') 
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Let  us  consider  (11)  and  substitute  x'^Ip)  from  (70)  for  the  range  0  to  ,  x(-’<|p)  ffom  (71)  for  the 
range-oo  to  0  and  x(x|33j)  =  Then  we  obtain 


C  i/)(x|p)  i^*(x'|p)  dp  +  r  ^('(xIp)  b*(p)  ^(^*(x'|-p)  dp 

do  do 

+  r  b(p)  i/-(x|-p)  i;>*(x'|p)  dp  +  C  |b(p)|^  iHx|-p)  iZ-^ix'l-p)  dp 

do  do 

pO  ^  ,  Wx|E.)  i^(x'lE.) 

+  \  it(p)|  ‘^'(xIp)  K’*(x'|p)  dp  +  I - 

-oo  i 


6(x  -  x') 


(72) 


We  note'at  once  that 


r  lb(p)|^  i/'(x|-p)  ^*(x'|-p)  dp  -  \  lb(p)|'^  !^(x|p)  i/d‘(x’|p)  dp 

Jn 


Also,  since  !b(p)|^  l  !t(p)|^  =  1,  we  can  add  the  fourth  and  fifth  integrals  of  (72)  and  obtain  as 
their  sum 

i-0 

\  ii'(x|p)  i|’’^(x'|p)  dp 

-OO 

Tills  integral  in  turn,  when  added  to  the  first  integral  of  (72),  gives  for  tlieir  sum 


^»00 

\  iKxlp)  i/'*(x'|p)  dp 
Finally,  wc  note  that  since  b(-p)  b'  (p), 

»oo  1*'* 

r  ii(x|p)  b*(p)  i^’"(x'|-p)  dp  +  \  b(p)  (/-(xl-j))  ;/'*(x’|p)  dp 

do  *'0 

.  r  ^(xlp)  b(-p)  li'^x’l-p)  dp  +  1°  b(-p)  ^'(xIp)  ^-"(xM-p)  dp 

do 


^»00 

\  il(x|p)  b(-p)  i/'’*(x'|-p)  dp 
«00 


With  these  simplifications,  (72)  now  becomes 

,.«o  P*  *  I  , 

\  if(x Ip)  (X' Ip)  dp  *-  \  i('(x1p)  b(-p)  ii<*(x'i-p)  dp 

d  .00 


^  nx\E.)  Wx'|E.) 

+  L  A. 

1 

1 


6{x  —  x') 


(73) 


It  now  appears  that,  as  a  consequence  of  (73),  we  can  expand  any  quadratically  integrable  function 
f(x)  in  the  following  manner; 


l'(x) 


i/'ixlp)  a(p)  dp  +  ^  i1)(x|E.)  a. 
i 


29 


where 


,P,  =  £ 

“I'i. 


fly)  !|'*(y|p)  dy  +  J 
Hy)  Hy\E^) 


f(y)  b(-p)  i/i^lyl-p)  dy 


A, 

1 


dy 


However,  it  follows  immediately  that 
a{p)  s  a(p) 

where  a(p)  is  given  by  (65).  Hence,  for  any  quadratically  intcgrable  function  f(x)  we  must  have 
for  all  p 


[b{-p)  -  b(-p))  r  f(x)  ^(i(x|p)  dx  E  0 

-CO 

A 

This  can  be  true  only  if  b(p)  =  b(p)  for  all  p. 

V.  SOME  EXAMPLES  OF  THE  SOLUTION  OF  INVERSE 
SCATTERING  PROBLEMS 


We  shall  now  give  some  examples  of  the  extreme  potency  of  the  Gel'fand-Lcvitan  algorithm 
and  show  liow  it  can  be  used  to  reconsti'uct  the  entire  scattering  problem  from  the  scattering 
coefficients. 


A.  The  Delta  Function  Potential 


First  let  us  consider  the  case  for  which  there  are  no  bound  states  and 


b(p) 


A  ® 
2 


-Ziop 


|3  +  ^ 


iA 

2 


(74) 


This  reflection  coefficient  is  of  the  proper  form,  since  it  satisfies  the  relations 
b(-|))  =  b*(p) 

lim  b(p)  =  0  in  the  upper  half  plane 

IpI  -«> 

In  this  case ,  we  now  have 


R(x) 


g-ip(2Q;+x) 


dp 


(75) 


Evaluating  the  integral  in  the  usual  fashion  by  contour  integration,  i.e.,  closing  the  contour  in 
the  upper  half  plane  when  x  <  —2a,  we  obtain 

R(x)  =  0  for  X  <  —2a 

Similarly,  closing  the  contour  in  the  lower  half  plane  when  x  >  —2a,  we  obtain 
R(x)=— -^e  when  x>— 2a 
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Hence, 


.  -Y(x+2a) 

R(x)  =  —  ^  i)(x  +  2a)  e 

The  Gel'fand-Lcvitan  equation  for  x  >  y  is  then,  for  the  given  example, 

A  -y(x+y+2a) 

K(x,y)  =  ^  77(x  +  y  +  2a)  e 

px  -^(z+y+2a) 

+  -y  \  K(x,  z)  r|(z  +  y  +  2a)  e  dz 


Wc  note  than  when 

X  <  —  a 
X  +  y  <  —  2a 
z  +  y  <  —  2a 

and 

I<(x,y)  =  0 

and  hence  also 


(76) 


(77) 


V(x)  0 

Equation  (77)  suggests  that  wc  write 

K(x,y)  =  i)(x  +  y  +  2a)  g(x,y)  .  (78) 

Substituting  tliis  expression  for  K(x,y),  wc  obtain 


-^(x+y+2a) 

T)(x  t  y  i  2a)  g(x,y)  =  t)(x  +  y  +  2a)  o 

px  -■^(y+zt2a) 

+ -—  \  T|(x  +  z  +  2a)  7)(y  +  z  +  2a)  g(x,z)  e  dz 


Consider  the  product  of  two  step  functions.  Evidently, 

[  T)(x  -  a)  when  a  >  b 


T)(x  -  a)  7)(x  -  b)  = 


I  i|(x  -  b)  when  b  >  a 

=  T)(a  -  b)  J|(x  -  a)  +  T)(b  -  a)  t)(x  -  b) 


Hence, 

T)(x  +  z  +  2a)  j;(y  +  z  +  2a)  =  T)(y  —  x)  r)(x  t  z  +  2a)  +  t)(x  —  y)  T)(y  +  z  +  2a) 

=  T)(y  +  z  +  2a) 

since  x  >  y. 

Using  this  expression  for  the  product  of  the  step  functions,  we  now  have 
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A  -y(x+y+2a) 

j)(x  +  y  +  2a)  g(x,y)  =  y  J}(x  +  y  +  2a)  e 


+ 


A 

2 


-y(y+z+2a) 

7)(y  +  z  +  2a)  g(x,  z)  e 


A  -y(x+y+2a) 

=  y  J)(x  +  y  +  2a)  e 


+  y  ))(x  +  y  +  2a)  J 


-(y+2a) 


-^(y  h:+Za) 
g(x,z)e  dz 


We  liavR  already  shown  that,  for  x  +  y  ^  2a  <  0,  K(x,y)  “  0.  Hence  we  now  treat  the  case 
X  +  y  +  2a  >  0  so  that  t)(x  +  y  +  2a)  =  1.  Then 


-y  (y+2a) 

=  o  f(x,y)  , 


whei'e 


f(x,y)  =  ^  o 


A  2^  ,  A 

+  y  e 


o'^^f(x,z)dz  , 


-(y+2a) 


and 


df(x,y)  _  A  ^A(y+a) 

-717-  -  T  ^ 


f(x,  -y  -  2a) 


We  ran  satisfy  (81)  by  letting 

A 

2^ 

f(x,y)  e  v(x) 


and  substituting  in  (80)  we  obtain 

A  A 

2^  ,  ,  A  ■  2^  A  -Aa  .. 
e  v(x)  =  y  G  7  ° 


A. 

2  '' 

e  dz 


-(y+2a) 


_  A 

A  2^  -Aa  ,  . 
=  y  c  -  e  v(x) 


-yX  y(y+2a)' 

o  -  e 


which  gives 


,  ,  A  Aa 
v(x)  =  y  c 


(79) 


(80) 


(81) 


Hence, 


f(x,y) 


A 

2  ® 
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and 


g(x,y) 


A 


Then 


and 


K(x,y)  =  7)(x  +  y  +  2a]  y 


V(x)  =  2  ^  [K(x,x)] 

=  A6(x  +  O') 

Having  determined  l•C(x•,y),  we  can  now  also  find  the  eigenfunctions  !ji(x|p): 

,ipx 


i/i(x|p)  = 


A 


Ztt  Z\l~2Tr  -<a 


r 


T)(y  +  X  +  Za)  e^'’^  dy 


.ipx 


A 


1)(X  +  Ci)  1 

•s/yir  2\f^  -{\+Za) 


e^*^^  dy 


ipx 


-  [1  -  7)(x  +  tt-)  I  +  J)(x  +  a)  e 


Ai 


^-2ipa  e_ 


•ipx 


Zp  2p 

r’rom  tliis  we  see  at  onee  that  i/'(x|p)  satisfies  (48),  namely, 

ipx 

lim  i/'(x|p)=-2- —  for  all  p 

X— -«  sfZTT 


\f2Tt 


Also, 


di/)(x||))  _  ip  e^‘^'''  ,, 
— - 


Ipx 


,  Ai  5(x  4  a)  e-2i>’“  " - 

4zi 


A  ,  ,  ,  -2ipQ-  e-'P^ 

+  7)(x  +  a)  e  - 


v  Ztt 


ip  o  ‘  '  ^  M  6(x  +  a)  ((,-2ipa  ^-ipx  _  ^ipx^ 
'/Ttt  ^P  \'~ZTr 

+  A  7i(x  +  a)  ^gipx  ^  g-2ipa  ^-ipXj 

2  4zi 


ip  e^P’^  ^  A  77(x  +  a)  ^^ipx  g-2ipa  ^-ipx^ 


'sIZn 


sfzH 


d^i))(x|p)  ^  _  2  ^  A  6(x  4-  a)  ,  ipx  ^  g-2ipa  ^-ipXj 


dx 


n/2^ 


'\mi 


ipA  i|(x  +  a)  ^gipx  _  g-2ipa  ^-ipXj 
2  n/27 


(82) 


(83) 


(84) 
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Hence, 


d^if'lxlp)  ^  _  2  e^P’^  ^  ipA  v(k  +  a)  ^^ipx  _  g-Zipa  g-ipXj  ^  g-ipa 


dx 


4zit 


\fZTi 


+V(x)  ^(x|p)  =  p^ 


dx 


.ipx 


Ai 


(1  -  ^  r)(x  +  a)] 

^/2^  ^P 


Ai  .  .  -Zina  e  '^P^  A  „-ipa 

+  7-i?(x  +  a)e  -—  -— =  e 

n/ZtT  I  's/Ttt 


ipx 


+  A6(x  +  o)  [1  -  ^  ^  '■  “> 


Ai 


-2ipa  e 


-ipx 


’■IzH 


P^i|'(x|p) 


i/i(x|p)  satisfies  the  Schrodinger  equation; 

i! 

2 


--^  +  A6(x  +  a)|  ijifxlp)  =  p^V'(x|p) 


dx 


(85) 


subject  to  the  boundary  condition 

gip’^ 

lim  i/ijxlp)  =  ——  for  all  p 

X-*-oo 

and  hence  for  p  >  0, 

X(x|p)  =  (/'(x|p)  +  b(p)  i/ifxl-p) 
is  tlie  unique  solition  which  satisfies  the  same  Schrbdinger  equation  subject  to  the  different  bound¬ 
ary  condition 


(86) 


lim  x(5^lp)  =  - -  + 


-ipx 


\;-*-oo  y  Zjt  Ztt 

Substituting  in  (86)  from  (74)  and  (84),  we  have 
X(x|p)  =  !p(xlp)  +  b(p)  i^(xl-p) 


ipx 


[1  -  ^  T)(x  +  a)]  +  ^  i)(x  t  a)  e 


Ai 


T  „"‘P^ 
-2ipa  e 

n/T^ 


^  e:^  M  _  A%(x  4  a)^ 

2  (p  +  lA,  VZ^  2p  4p(p  +  i  f) 


ipx 


Vzlr 

.  e-iP’^  I  Ai 


,  Aiii(x  +  g)  _  A  n{ 

2p 


4p(p 


X  +  g )  I 

X  •  A, 

4  i  T )  J 


Ai  ^  2d  ^  -Ziap 

^V(x  +  a)-^^  - « 

(p  +  -^) 


^/2¥  pP 

Letting  x  —  "  ,  we  can  now  also  find  the  transmission  coefficient  t(p) 
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lim  x(x1p) 

X-»oo 


1  _  Ai  A^ 

Ai  Ai  2p' 

4^ 

^  4p(p  +  if^)J 

sfzir 

[2p  2p 

^-Ziap 


pipx 

lim  x(xlp)  =  t(p)  — - 

x-*«>  nZH 


where 


t(p)  =  1 


M 

2p 


2p(2p  +  iA) 


or 


t(p)  = 


_^12_ 

2p  +  lA 


(87) 


Also,  as  has  been  shown  to  be  generally  true, 

""  -2p  +  iA  ^  2'P^TA  " 

and 

|b(p)|^  +  |t(p)|^  2+  ■ 

+  A  4p  +  A 

We  can  also  find  r(p),  the  reflection  coefficient,  from  the  other  side  by  proeeeding  in  the  following 
manner. 

For  p  <  0, 


lim 

X->-oo 


X(x1p)  =  s(p) 


n/27- 


It  was  shown  earlier  in  this  report  that 

s(p)  =  t(p)  for  all  p 


Hence,  we  write  for  p  <  0; 

X(x1p)  =  s(p)  i|)(x1p) 


2P 

2p  -  iA 


[  1  1 


+  a)  e 


-2ipQ 


42^ 


In  order  to  find  r(p),  we  now  let  x  +°°  and  find 


lim  x(x|p)  = 

x-*«> 


^P, 

2p  —  iA 


,1  _  M  e-^'P“ 

'  2p'  2p  ^ 


e’P^ 


e^P^  Ai  e“^^P“  ^ 

^  2p-Ai  ^ 


35 


Hence, 


r(p) 


Ai  -Zipoi 
2p  -  Ai 


We  might  note  that  only  for  a  =  0  do  we  have  r(p)  =  b(p)  for  all  p. 

Before  proceeding  to  further  examples,  we  emphasize  again  the  extreme  potency  of  the 
Gel'fand-Levitan  algorithm.  In  the  previous  example,  we  have  reconstructed  the  entire  scattering 
problem  from  a  knowledge  of  the  reflection  coefficient  on  one  side.  We  have  found  the  potential 
function  producing  this  reflection  coefficient,  the  corresponding  eigenfunctions  of  the  SchrSdmger 
equation,  the  transmission  coefficient,  and  the  reflection  coefficient  on  the  other  side. 


B.  Reflectionless  Potential 

In  this  example,  we  shall  \ise  the  Gel'fand-Levitan  algorithm  to  construct  a  reflectionless 
potential. 

Wc  put  b(i))  i  0,  and  let  us  also  assume  that  there  is  but  one  bound  state,  L^,  having  a  noi-- 


malizat  ion  A  j .  'I'hen 


R(x) 


(88) 


n'  ‘"1 


(x+y) 


K(x,y)  =  -  ■ 


A. 


-i': 


.y-Hi(z+y) 


K(x,  z) 


dz 


K(x,  z)  e 


dz 


(89) 


Equation  (89)  suggests  tliat  wc  write 


K(x,y) 


A, 


f(x) 


(90) 


where 


-E.x 


f(x)  e' 


fLo  r 
Ai  J. 


X  2j-E,z 


dz 


f(x)  e 


(91) 


(92) 


and 
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(93) 


K(x.y) 


+  2AJ-Ejc 


-Ejx 


••  rr-  I — 

=  ®  soch  J-E^  (x  -  x^)  , 


where 


In  2A,J-E, 


Hence,  \vc  obtain  a  reflectionloss  pott'nlial  \'(x)  given  liy 


Also, 


i;.(x|i))  -  -5^ 


^  ^  c'fe  soch^y^Ej  (x  -  x^) 

_,ipx  l-\r^  ox  /-^(y-xj 


n/2^ 


y 


Serb  (X  -  .x^)  e'*’-''  dy 


gipx  _ _  .soch  /-E  (X  -  X  )  c'''  ‘  “ 

— -^'^1 - - 


S2- 


,ipx 


*■  ip) 


1  - 


v/-l-l  =^ech^-E,  (x-x^^)  o 


Since  b(p)  0, 


e^'-- 

vTrr 


^/-E,  scchJ-E,  (x-x^^)  c 


J-E,(x-x^)-| 


/-'■'l  '•  ip 


In  order  to  dolerniine  t(p),  we  ron.sidor 


II  once. 


lim  x(xip) 

\  ►  CO 


t(p) 


e‘PX 


2.  /-  !■: . 


'-ip 


p-iJ-E 


We  may  now  write  for  p  <  0, 

,ipx  (p  -  i^-E^) 


X(x|p) 


xTzi  (p  +  iJ-E^) 


1  - 


/-E^  sechJ-Ej  (X  -  x^)  o 


/  i^i(x  x^)^ 


-Ej  +  ip 


(9d) 


(95) 


(96) 


(97) 


(98) 


(99) 
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Then 


lim 

X-^M 


X(x1p)  — 

42t! 


and  hence  r(p),  the  reflection  coefficient  on  the  other  side,  is  also  zero.  Tliat  this  must  be  so 
is  apparent  at  once  from  the  relations 


Finally , 


lb(p)l^  +  |t(p)l^  =1  for  p  >  0 

|r(p)|^  +  |s(p)|^  =1  for  p  <  0 

s(p)  =  t(p)  for  all  p 
wo  might  mention  that  in  this  case 


X(x|  E^)  =  i()(x|E^) 


sech 


(x  -  x^)  dv 


sechJ-E^(x-x^) 


/- 


E,(2x-x^) 


sech^^E^  (X  -  x^) 


(100) 


Again  the  entire  scattering  problem  has  been  reconstructed  witli  the  help  of  the  Gcl'fand-Lcvitan 
algorithm. 


VI.  REFLECTION  COEFFICIENTS  FOR  POTENTIALS 
WHICH  VANISH  IDENTICALLY  FOR  x  <  -a 

Unfortunately,  the  potentials  we  have  been  obtaining  so  far  are  not  of  practical  interest,  since 
they  do  extend  to  infinity  and  are  therefore  impossible  to  construct  in  actuality. 

Wo  shall  next  give  conditions  on  the  reflection  coefficient  b(p)  which  will  result 'in  potentials 
which  will  be  identically  zero  on  one  side  at  least. 

Wo  first  prove  tlie  following  theorem; 


Thoo"em 


It  U(x)  a  0  forx<-2a 


Then  V(x)  ==  0  forx<-a 

Wo  hope  to  be  able  to  clioosc  b(p)  and  the  point  eigenvalues  E.  so  tliat  ll(x)  =  0  wlien  x  < -2n  and 
thus  obtain  a  corresponding  potential  V(x)  which  is  identically  zero  for  x  <  -  a.  To  prove  the 
theorem,  we  use  the  G el 'fand- Levitan  equation  for  x  >  y: 


K(x,y)  -R(x  +  y) 


K(x,  z)  H(y  +  z)  dz 
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If 


then 

and 

hence. 

Also, 

and  hence 

and 

Therefore, 

and 


X  <  —  O'  , 

y  <  -a 

X  +  y  <  — 2a 

R(x  +  y)  i  0 

z  <  X  <  —  a  , 
y  t  z  <  —  2n' 

K(z  +  y)  :  0 

K(x,y)  =  0  when  x  <  —  a 

V(x)  ^  2  ^  [K:(x,x)]  =  0  when  x  <  -  a 


This  eompletes  the  proof  of  the  theorem. 

Let  us  now  consider  b(p)  in  the  complex  plane.  Suppose  b(p)  has  polos  at  p  =  ITj  where  Tj  >  0. 
We  shall  show  that,  if  Kj  =  and  Aj  is  chosen  properly,  we  will  be  able  to  make  R(x)  =  0  for 
X  <  -2a;  hence,  b(p)  will  be  the  rofleqtion  eoeffieiont  for  a  potential  V(x)  which  is  identically  zero 
for  X  <  -  a . 

Assume  that  b(p)  has  the  form 


b(p'  -  g(p)  e‘^‘“'’  , 

where 

g(-p)  =  B*(P)  on  ncftl  axis 

-  k*  4-1 

g(p)  -  0(p  )  whore  K  >  1 

and  g(p)  has  poles  of  residue  r^  at  the  points 

p  -  ).r  =  i  Pe.  >  0 
‘  J  J 


(101) 


(102) 


(103) 


Then 


R(x) 


g(p)  dp  .  E 

j 


(104) 
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Wlien  X  <  -2a, 


1  r  ,  ,  -ip(x+2a)  ^  .  V 

^  \  g(p)  e  '  dp  =  1  2,  1'  e  J  J 

-00  J 


J 


V  ■  3  1 

2  ir.  e  •'  e  ■' 


2aT.  /-E.x 


Idenco,  if  we  choose  A.  so  that 
3 


2aT. 


A.  =  -'‘’j  ° 

J  ■’ 

tlnen  R(x)  =  0  for  x  <  —  2a,  and  hence  V(x)  s  0  for  x  <  —a. 
As  an  example,  consider  the  case  wliorc 


b(i))  = 


iB  e 


-2iap 


2  iB 

P  -  T 


In  this  case. 


77  \  >’(P)  ‘■'P  =  ^  j 


->1/3  ^-ip(x-t'2a) 


iB 

p--r 


dp 


Hence, 


forx<-2a 


K(x)  -  -il  c'^“  o“  r  — , —  for  x 

M 


-2a 


If  we  clioose 


1  B  B(v 

~  2  ° 

•M 


(105) 


(106) 


(107) 


(108) 


tlien 


K(x)  E  0  for  X  <  -2a 

Wlicn  X  —2a’,  we  close  the  contour  used  to  evaluate  tlic  integral  in  tlic  cxpi’cssion  for  ll(x)  in 
tlie  lowoi’  half  plane.  This  contour  contains  no  |)olos  of  tlio  integrand,  and  hence  in  tliis  case  tlie 
value  of  Hie  integral  is  ■/'oro.  Tlien 


R(x) 


A., 


B 

2 


B 

2 


(x+2a) 


for  X  >  —  2a 


(109) 
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Wc  now  have 


„  f-(x+2ar) 

R(x)  =  £  7}{x  I  2a)  e 


and  for  x  >  y. 


R  §(x+y+2a) 

K(x,y)  =  -  ^  7](x  +  y  +  2a)  e 


-if 


F<(x,z)  t)(y  +  z  F  2a)  c 


2  (y+z+2a) 


(Iz 


wlicn 


and 


and  hcnco. 


X  <  —  a 
X  F  y  <-  —  2  a 
z  r  y  <  —Za 

K'x.y)  .  0  , 


V(x)  =  0  . 

Equation  (HI)  suggests  that  we  write 

K(x,y)  =  i!(x  +y  t  2a)  g(x,y) 

and,  proceeding  as  in  a  previf)us  example  given  in  tliis  report,  wo  obtain 

fi(x,y)  =  -^^  , 


K(x,y)  =  -  -2  ’((X  '  y  >■  2a) 

V(x)  ^  -H6(x  +  o) 

As  wo  expected,  tliis  potential  is  identically  zero  for  x  <  —  q',  Suppose  now  that 
,  ,  ,  iB  1 

2 - m 

P-- 

and  that  there  are  no  point  eigenvalues.  Wo  will  now  not  be  able  to  make  ri(x)  s  0  for  x 
Let  us  see  what  the  corresponding  potential  will  be  in  this  ease. 

We  have 

H  , 

U(x)  -1*^^  ■ 

.0  for  X  >  0 


R(x)  =  -  —  T](-x)  e  ^ 


(110) 

(111) 


(H2) 

(113) 

(114) 
(H5) 

(116) 

<  — 2a . 


(117) 
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Hence,  for  x  >  y, 


K(x,y)  =  ^  7j(-x  -  y)  e  f  J  K(x,  z)  (-y  -  z) 

The  solution  of  this  equation  may  be  shown  to  be 


|(y+z) 


dz 


K(x,y) 


y(x+y) 

Byj-x)  e'^  ,  B 


2  —  e 


Bx 


+  ^  7](x)  7)(x  +  y)  +  Briix)  r){~x  -  y)  e 


B/  ,  ^ 
T(^+y) 


Then 


and 


K(x,x)  =  Bq(-x)  e 
2  -  e 


Bx 


t  -0  7)(x) 


V(x)  =  2  ^  [K(x,x)]  4B^.|(  X)  e^^^  ^ 


(2-0^^ 

^  4B^t)(-x)  _ 

(2-0^^^ 


B6(x) 


Hence,  in  this  case  the  potential  has  an  exponential  tail  for  x  <  0. 
Finally,  let  us  consider  as  an  example  the  case  for  which 


b(p) 


(p  4  i)  (p  -  i) 


E,  =-l 


In  this  case , 


-ip(x+2«)  X 

iMiid.--)  *  *  a; 


(118) 


(119) 


(120) 


(121) 


(122) 


(123) 


When  X  <  -2a,  the  exponent  in  the  integral  is  positive,  and  we  close  the  contour  in  the  upper  half 
plane.  Evaluating  the  integral  thus,  we  obtain  for  x  -2o: 


^  pw  -ip(x+2a) 


.J_  P  _£ 

’  27r  J  „  (p  + 


(p  1-  i)  (p  -  i) 


dp  =  —  i 


^x+2a 

2i 


-  e 


X  e 


2  a 


If  wo  choose 


1  c 


A 


(124) 


then  R(x)  0  for  x  <  —  2a,  and  from  the  theorem  wo  expect  V(x)  =  0  for  x  <  —  a.  When  x  —2a, 
the  exponent  in  the  integral  is  negative,  and  wo  close  the  contour  in  the  lower  hall  plane.  Thus 
we  obtain  for  x  >  —  2a; 


1  f"  -ip(x+2a) 


-si 


27ri  e 


-x-2a 


-2a 


(p  +  i)  (p  -  i)  (-2i) 


•  =  — e 
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Hence, 


or 


F’or  X  >  y , 


-2a  ,  2a 

R(x)  =  -e“^  ^ =  sinh(x  +  2a)  for  x  <-2a 


H(x)  =  ?)(x  +  2a)  sinh(x  +  2a) 


K(x,y)  =  -i)(x  +  y  +  2a)  sinh(x  +  y  +  2a) 

»x 

K(x,z)  7)(y  +  z  +  2a)  sir.’i(y  +  z  t  2a)  dz 

-OO 


We  note  that  when 

X  <  —  a 
X  +  y  <  —  2a 


z  +  y  <  —2a 


lienee 

K(x,y)  =  0 

and 


V(x)  =  0  when  x  <  —  a 

as  we  expected.  Let  us  now  try  to  solve  tlie  integral  equation  (126)  and  find  V(x)  foi-  x  >  - 
we  write 


K(x,y)  =  i)(x  +  y  +  2a)  g(x,y) 


tlien 

j)(x  +  y  +  2a)  g(x,y)  =  -?)(x  +  y  +  2q)  sinh(x  +  y  +  2a) 


g(x,z)  i)(x  +  z  +  2a)  J)(y  i  z  +  2q)  sinh(y  +  z  +  2a)  dz 


Hence, 


=  -T7(x  +  y  +  2a)  sinh  (x  +  y  +  2a) 


—  r  g(x,  z)  i)(y  +  z  +  2a)  sinh(y  +  z  +  2a)  dz 

-OO 

-T7(x  +  y  +  2a)  sinh(x  +  y  +  2a) 


-i?(x  +  y  +  2a) 


r 


(y+2a) 


g(x,z)  sinh(y  +  z  +  2a)  dz 


g(x,  y)  =  —  sinh  (x  +  y  +  2a)  —  ^  g(x,  z)  sinh  (y  +  z  +  2a)  dz 

‘^-(y+2a) 


(125) 

(126) 


a .  It 

(127) 


(128) 
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Also, 


t 


=  —cosh  (x  +  y  +  2a)  —  T  g(x,  z)  cosh  (x  +  y  +  2a)  dz 

'^-(y+2a) 


d  g(x,y) 
dy 


=  -sinh  (x  +  y  +  2a)  -  1  g(x,  z)  sinh(x  +  y  +  2a)  dz  -  g(x,  -  y  -  2a) 

2  '^-(y+2a)  ' 


or 


=  g(x,y)-g(x,  -y-2a) 

dy^ 

This  is  satisfied  by 

g(x,y)  =  v^(x)  u(y)  +  v^lx) 

where  ii(y)  satisfied 

d^u(y) 


dy' 


From  (131). 


and 


u(y)  -  u(-y  -  2a) 

u(y)  -  sinh  sTz  (y  t  a) 


(129) 


(130) 


(131) 


(132) 


g(x,y)  =  v^(x)  sinh '\/^(y  +  a)  t  v^(x)  .  •  (133) 

Substituting  this  expression  for  g(x,y)  in  the  integral  equation  (128),  wo  obtain 
v^(x)  sinhx^  (y  +  a)  +  v^ix)  =  -sinh  (x  +  y  +  2a) 

-  C  (v,(x)  sinh\f2  (z  4  a)  +  v..(x)]  [sinh(y  +  z  +  2n)l  dz 

d-(y+2a)  ^ 


=  -sinh(x  !  y  +  2a)  -  v^ix)  lcosh(y  +  z  (■  2a)l^jy^2n') 


-v^(x) 


2(1  +  n/2) 


-  sinh  [(1  +  \f2)  z  4-  y  +  (2  +  V  2)  a  ] 


- - -  sinh  [(1  -  x/T)  z  4- y  4-  (2 xf2)  a  1 1 

2(1  -V2)  .  I  -(yt2a) 

=  -sinh(x  4-  y  4-  2a)  -  V2(x)[cosh(y  4-  x  4-  2a)  -  1| 

1 


-v^(x) 


2(1  4-  \/2) 


sinh  [x  4-  y  4-  2a  4-  (x  4-  a)] 


j 

+ 

sinh^/?  (y  4-  a) 

2(1  4-  \f2) 

1 

sinh  [x  4-  y  4-  2a  —  (x  4-  a)] 

2(1  -  4z) 

1 

4- 

sinh^/?  (y  4-  a)j 

2(1  - 

1 
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or 


0  =  —  sinh  (x  f-  y  +  2a)  —  v^lx)  cosh  (x  +  y  +  2a) 


_V  (X)  - — 

1  2(1  +  ^fz) 


Sinh[(x  +  y  +  2q)  +  sfz  (x  +  a)] 


2(1  -  V2) 


sinh  [x  +  y  +  2q  —  s/Y  (x  +  a)] 


0  ^  _^x+y+2n  ^  ^-(x  ly.+2cv)  _  j^x+y+2«  ^  p-(x+y+2cv)j 


j^x+y+2Q+'/2(xiQ)  _  g-(x+y+2Q)-"/2(x+0')j 

^  2(1  +  n/2) 


1 


2(1  -^|Z) 

^x+y+2  Q' 


,  x+y+Z  a  -4z{x+a)  _  -{x+y+Za)+'4z{x+a), 

1 G  G  J 


,  ,  \IZ{x+a)  I  ,  -\IZ{x+a)  , 

v,(x)  e  '  V  (x)  e  '  ^ 

-1  -  v,(x) - - -  +  — - 

2(1  +  n/7)  2(1  -  \[Z) 


+  e 


-(x+y+2Qi) 


1  -  v.{x)  + 


v.j(x)  e 


-4z{x+a) 


v^(x)  e 


\^(x+q) 


2(1  +  \fZ) 


2(1  -  s/2) 


Hence,  Vj(x)  and  V2(x)  must  satisfy 


1  +  V2(x)  + 
1  -  v^ix)  - 


v.j(x) 

^s/2(x+cv)  ^-s/2(x+a) 

=  0 

2 

i  1  +  s/2 

1  _  vl 

Vj(X) 

■  ^42(x+a)  ^-s/2(x+qi) 

=  0 

2 

l-s/2 

1  +  s/2 

Therefore, 

^1*-'^^  1  ^^(x+o;)  [__! _ 1 _ )  ^  3-n^(x+«)  [_1 _ 1 _ 

^  2  I®  U  +  s/2  l-sT^J  h  +  s/I  1  -  s/2 

.2+v^(x)-/2[e^<-'^+“>+e-^<^+“>l  =  0  , 


or 


Also, 


and 


V^(x)  =  - 


1 

42 


_ 2 _ 

j^s/I/x+it),  ^-42{xi-a) 


_ 1 _ 

4z  eoshs/2(x  +  a) 


^^4z{x+a)  _  ^-4z{x+a)^_^^^4Z{x+a)  ^  ^-4z{x+a)^  _  ^ 

1  +  Vt(x)  +  ! - — - - - —  ~ 

s/2 


(134) 


(135) 


& 
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(136) 

(137) 

(138) 

(139) 

(140) 
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